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Superstatistics are superpositions of different statistics relevant for driven nonequilibrium systems 
with spatiotemporal inhomogeneities of an intensive variable (e.g., the inverse temperature). They 
contain Tsallis statistics as a special case. We develop here a technique that allows us to analyze the 
large energy asymptotics of the stationary distributions of general superstatistics. A saddle-point 
approximation is developed which relates this problem to a variational principle. Several examples 
are worked out in detail. 
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I. INTRODUCTION 

Many driven nonequilibrium systems exhibit complex 
dynamical behavior characterized by spatiotemporal fluc- 
tuations of an intensive parameter [3 which may represent 
an inverse temperature, a chemical potential (e.g., in a 
system with inhomogeneous concentrations), an effective 
friction constant, the amplitude of a perturbing noise, or 
the local energy dissipation, as in the case of turbulent 
flows. Although such systems do not settle to equilib- 
rium, their long-term behavior can often be described in 
the spirit of equilibrium statistical mechanics by view- 
ing them as consisting of an ensemble of subsystems or 
"cells" to which are associated different values of (3. If 
there is local equilibrium in each cell, so that statisti- 
cal mechanics can be applied locally, and if the fluctua- 
tions of /3 evolve on a sufficiently large time scale, then 
in the long-term run the entire system can be described 
by a mixture or superposition of different Boltzmann fac- 
tors having different values of (3. Such a mixture of vari- 
ous statistics has been termed a "superstatistics" and 
has been the subject of various papers lately (see, e.g., 
Refs. 0, S i, m II II S in El d) Many mod- 
els based on the notion of superstatistics have also been 
applied successfully to a variety of physical problems, 
including Lagrangian 0, [l3 and Eulerian turbulence 
[isl ll^ , defect turbulence cosmic ray statistics ^3 j 
plasmas statistics of wind velocity differences |20j . 
and econophysics 0, j^^ . 

What is common to all these problems is the exper- 
imental observation of stationary distributions having 
"fat" tails. Such distributions fall necessary outside the 
framework of ordinary statistical mechanics, but not that 
of superstatistics. For example, if the random intensive 
parameter /3 in the various cells is taken to be distributed 
according to a particular probability distribution, the x^- 
distribution, then the corresponding superstatistics, ob- 
tained by integrating the Boltzmann factor e~^^ over all 
/?, yields the nonextensive statistics of Tsallis defined by 
the so-called g-exponential function [231 123. 125ll2^ 

e-^"^-[l + (g-l)/?oi?]-'/(^-'^ (1) 

This particular statistics decays as a power-law for large 
energies E rather than an exponential, as is the case for 



the ordinary Boltzmann factor. In this sense, it is a fat- 
tailed statistics. The parameter /3o above is related in the 
superstatistical model to the average inverse temperature 
of the inhomogeneous system, whereas the so-called cn- 
trop ic index q relates to the variance of the (3 fiuctuations 
[23,123. It is worth mentioning that distributions having 
the form of a g-exponential can be obtained formally by 
maximizing Tsallis' measure of entropy subject to suit- 
able constraints. Moreover, ordinary statistical mechan- 
ics, which correspond in the superstatistics picture to the 
case where there is no fiuctuations in (3, is recovered in 
the limit where g ^ 1 [2l0,0,|2|. 

For other distributions of the intensive parameter /3, 
one ends up with more general superstatistics. Gen- 
eralized entropies, which are analogs of the Tsallis en- 
tropies, can also be defined for these general superstatis- 
tics 1^ ^3 J ^'^d generalized versions of statistical me- 
chanics can be constructed, at least in principle. It has 
been shown that the corresponding generalized entropies 
are stable |llTl|. 

In this paper we will briefly review the superstatis- 
tics concept, and then analyze the asymptotic behavior 
of general superstatistics for large values of the energy 
E. We will investigate how the properties of the func- 
tion /(/3), which represents the probability distribution 
of the intensive variable (3 in the various spatial cells, 
determine the asymptotic decay rate of the generalized 
Boltzmann factor of the superstatistics. We develop a 
saddle-point approximation technique which allows us to 
treat this problem in full generality. Several examples 
will be worked out in detail to show that the asymptotic 
decay rate of the stationary distributions resulting from 
f{(3) can not only be a power-law, as in the case of nonex- 
tensive statistics, but can also be an exponential of the 
square root of the energy or, generally, a stretched ex- 
ponential. We will discuss universal aspects of the large 
energy asymptotics, thus complementing the considera- 
tion in Ref. |lj] where universal aspects of the low-energy 
behavior of general superstatistics were discussed. The 
large energy asymptotics is of particular physical impor- 
tance because the tails of the observed distributions mea- 
sured in various experiments (e.g., hydrodynamic turbu- 
lence '2^, plasmas IstI, and granular media Sl]) can 
distinguish between the various possible types of super- 
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statistics. 



II. SUPERSTATISTICS: BASIC CONCEPT 

Let us first briefly review the superstatistics concept 
as introduced in Ref. (Ij. Consider a driven nonequi- 
librium system witfi spatiotemporal fluctuations of an 
intensive parameter /3, for example, the inverse temper- 
ature. Locally, i.e., in spatial regions (cells) where /3 is 
approximately constant, the system is described by or- 
dinary statistical mechanics, i.e., by an ordinary Boltz- 
mann factor e~^^ , where E is an effective energy in each 
cell. In the long-term run, the system is described by 
a spatiotemporal average over the fluctuating /3. In this 
way, one may define an effective Boltzmann factor B{E) 
for the whole system as 



1*00 

B{E) = / /(/3) 



(2) 



where /(/?) is the normalized probability distribution de- 
scribing the /3- fluctuations in the various cells. For so- 
called type-A superstatistics, one normalizes this effec- 
tive Boltzmann factor and obtains the stationary, long- 
term probability distribution. 



where 



p{E) - ^B{E), 



Z= I B{E)dE. 
'0 



(3) 



(4) 



For type-B superstatistics, the /3-dependent normaliza- 
tion constant of each local Boltzmann factor is included 
into the averaging process. In this case, the invariant 
long-term distribution is given by 



00 ^-pE 



d/3. 



(5) 



where Z{j3) is the normalization constant of the Boltz- 
mann factor e~^^ for a given /3. Both approaches can 
be easily mapped into each other by considering a new 
probability density /(/3) = c- f{/3) /Z{P), where c is a nor- 
malization constant. One immediately recognizes that 
type-B superstatistics with / is equivalent to type-A su- 
perstatistics with /. 

As mentioned before, the fluctuations of the intensive 
parameter (3 are spatiotemporal: they can be produced 
by either temporal changes of the environment or by the 
movement of a test particle through inhomogeneous spa- 
tial regions. The precise nature and behavior of P in these 
situations can be quite varied. In our general description 
of superstatistics, we have taken P to be an inverse tem- 
perature which varies randomly in time or in space, but 
(3 can also represent, say, a chemical potential that varies 
smoothly in space. In this way, one may study supersta- 
tistical models where the fluctuations of P are caused by 



large-scale temperature gradients in a system [3 or by a 
nonuniform chemical potential describing inhomogeneous 
concentrations spread in space. 

From a more dynamical perspective, a superstatistics 
can also be achieved by considering Langevin equations 
whose parameters fluctuate on a relatively large time 
scale (the adiabatic regime [s^l); see Ref. 0| for details. 
For turbulence applications, for example, one may con- 
sider a superstatistical extension of the Sawford model 
of Lagrangian turbulence 0, 0, [s^l • This model con- 
sists of suitable stochastic differential equations for the 
position, velocity, and acceleration of a Lagrangian test 
particle in the turbulent flow. The superstatistical ex- 
tension naturally enters due to the fact that the local 
energy dissipation rate in turbulent flows is a random 
variable. Thus the parameters of the Sawford model be- 
come random variables as well. These kinds of models 
well reproduce experimentally measured turbulence data 
^13, 14]. 

In the following we will use the notation of type-A 
superstatistics, keeping in mind that we can always pro- 
ceed to type-B superstatistics by replacing / by /. We 
will restrict ourselves to positive values of P and E and 
will assume, in addition, that f{P) is everywhere differ- 
entiable and unimodal (single-bell-shaped curve). 



III. LOW ENERGY ASYMPTOTICS 

We recall the low-energy asymptotics of superstatistics 
for reasons of completeness; it was previously discussed 
in Ref. P]. Consider a distribution /(/?) having mean 
(P) — Po and variance 

(P^) - iPf = (P^) - P^, = a\ (6) 

Using the definition of B{E), we can write 

BiE) = {e-P^) = e-'3oS (^e-(P-P«)E^ . (7) 

Then, expanding in Taylor series the exponential term 
inside the expected value, we obtain 

B{E) = e-^»^ f 1 + \a^E^ + {iPo - Pf) ^\ • 



/c=3 



(8) 



Thus, to second order in E, B{E) must behave like 

B{E) ^ e-'^"^ ( 1 + -<j^eA (9) 



as E —t 0. This approximation represents the leading 
order correction to ordinary statistical mechanics in our 
nonhomogeneous system with temperature fluctuations 
for small values of the energy E. The zeroth-order ap- 
proximation to B{E) corresponds, as is expected, to the 
"pure" Boltzmann statistics. 



B{E) - e 



-l3oE 



(10) 
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with inverse temperature (3q. It can be noted that 

these two asymptotic results can also be considered to fiP) 

be valid approximations of B{E) in the limit where 

((/3o — for all fc > 2, that is essentially when 

/(/3) — > S{l3 — Po) (small fluctuations limit). 



IV. HIGH ENERGY ASYMPTOTICS 

To find the high-energy asymptotics of B{E), we use 
the fact that the integral defining B{E) has the form of a 
Laplace integral for £' — > cx) [2J|. In this limit, the inte- 
gral can be approximated by its largest integrand. This is 
the essence of Laplace method, otherwise known as the 
saddle-point approximation 0) ■ The conditions of 
applicability of this approximation method are basically 
the conditions that we have assumed regarding the shape 
of /(/?) and its differentiability. Thus, putting B{E) in 
the form 



B{E) = 



-/3S+ln/(/3) 



d(3, 



(11) 



we attempt to locate the largest integrand by finding the 
unique value of /3 which maximizes the exponent func- 
tion. 



<i>{(3,E) = -f3E + \nf{/3), 



(12) 



for any large enough energy value E. We call the value of 
f3 maximizing $(/?, E) for fixed E the dominant inverse 
temperature and denote it by (3^- The fact that /(/?) is 
assumed to be unimodal ensures us that (3e is unique, as 
required. Indeed, observe that ln/(/3) must be a concave 
function of (3 if /(/3) is unimodal (see Fig. 1), and, in this 
case, the maximum of $(/3, E) along the /3-direction can 
only be attained at a single point Pe which is such that 



i?- [In /(/?)]' 



f'iP) 



(13) 



Solving this equation for we find Pe and thus write 



BiE) 



^-PsE+lnfifSE) 



fWE)e- 



(14) 



in the limit where E ^ oo. Note that this basic Laplace 
or saddle-point approximation of B{E) can be improved 
a bit more by evaluating the integral defining B(E) using 
a Gaussian approximation of the integrand 34]. What 
results from this more refined approximation is the fol- 
lowing high-energy asymptotics: 



B{E) 



f{pE)e- 



-I3eE 



V-[ln/(/3B)]" 



(15) 



which differs from the Laplace approximation only by the 
square-root term involving the curvature of ln/(/3). 



ln/(/3) 




FIG. 1: Plot of a typical unimodal distribution /(/3) (top) and 
its logarithm (bottom). The dominant inverse temperature 
I3e is such that In /(/3b) has a slope equal to the energy value 
E. 



The Laplace approximation of B{E) as well as its 
Gaussian corrected version shown above are quite inter- 
esting from the physical point of view because they show 
that the mixture of Boltzmann statistics defining B{E) 
reduces at high energy to a "pure" Boltzmann statis- 
tics, just like in the equilibrium situation, although now 
the Boltzmann statistics involves an inverse temperature 
(3e which depends on the energy E considered. This 
means that, for high values of E, the long-term, station- 
ary behavior of the noncquilibrium system considered is 
dominated by the equilibrium behavior of a subset of 
cells having an inverse temperature equal or close to (Se- 
How Pe changes as a function of E is determined from 
the properties of /(/3). In the general case considered 
here, where ln/(/3) concave, we find that Pe when 
E oo (see Fig. 1). Thus in this case the large-energy 
behavior of the generalized Boltzmann factors B{E) is 
determined by the small-/? (high temperature) behavior 
of the function f{P). 

To complete this section, note that the exponent func- 
tion ^{Pe,E) which enters in the asymptotics of B{E) 
represents nothing but the Legendre transform of In /(/3). 
The result of this transform is a function of E which can 
be thought of to represent, following the theory of large 
deviations [3^ . an entropy function if we consider that 
the function ln/(/3) represents a free energy function. 
This entropy function, however, is just a formal one and 
is unrelated to the Tsallis en trop ies or other generalized 
entropies as defined in Refs. pXl l23|. 
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EXAMPLES 



B. Exponential tail 



We now consider a few cases oi (3 ^ asymptotic be- 
havior of f{P) and derive the corresponding asymptotic 
behavior of B{E) ior E ^ oo using the Laplace approxi- 
mation of B{E) corrected with the square-root term, i.e., 
Eq.lUni). 



A. Power-law tail 

Consider first an /(/3) with ^ , 7 > for 

/3 — > 0. An example of probability density having this 
asymptotic form is the following distribution for (3 



/(/3) = 



1 



r(§) V2/3o 



n/2 



(16) 



iPo > 0, n > 1), which behaves as /(/?) ~ /3"/^ ^ around 
/3 = 0. Another example is the F-distribution [llIT^. 



fiP) = c 



(1 -|_ 2^p){v+w)/2 ' 



where v,w,b are parameters and C is a normalization 
constant. 

To find the asymptotic form of the superstatistics 
B{E) corresponding to the choice ~ fi'^ , we pro- 

ceed to determine the value of the dominant inverse tem- 
perature (3e by solving Ea. lfT^ . The solution here is 

/3_E ~ l/E, so that 



l3EE + \nf[pE) 



-7 + 7 In 7 — 7 In 
--l\nE 



and 



[ln/(/3ij)]" 



-E' 



Combining the two results in H14|l . we obtain 



B{E) 



E 



= E-^-\ 



(18) 



(19) 



(20) 



Thus we see that power laws in /3 for small /3 imply a 
power law in E for large E, no matter what the rest of 
the distribution /(/3) looks like. Comparing this asymp- 
totic result with the g-exponential distributions studied 
in nonextensive statistical mechanics |3 HE HEli 
which asymptotically decay as -^^g ggg that 



7 + 1 



1 



(21) 



Power-law superstatistics are physically relevant for 
many different physical problems: e.g., defect turbulence 
cosmic ray statistics ^^d wind velocity mea- 
surements j23| among others. 



Consider now a density /(/3) having the asymptotic 
form /(/3) ^ e.~'^^^ , c > 0, as /3 — > 0. An example is the 
inverse distribution 



/(/3) = 



/3o 

r(f) 



V2-2. 



(22) 



(/3o > 0, 71 > 1), which behaves, for ^ 0, as follows: 



e 2/3 



(23) 



This form of /(/3) was previously considered in the con- 
text of density fluctuations in fusion plasma experiments 
M| as well as temperature fluctuations in perfect gases 
and arises if the temperature T = {kBP)~^ rather 
than P itself is x^-distributed. 

For this example, the equation that we need to solve 
to find /3b is simply 



[ln/(/3)r = ^ 



E, 



(17) and so (3e — 



'c/E. As a result, 

-PEE + \nfiPE) = -2V^ 



and 



BiE) 



-0EE+\nf{l3E) 



V-[ln/(/3B)]' 
using the fact that 

Pe 



2^3/2 



-2VcE 



-E^'\ 



(24) 



(25) 



(26) 



(27) 



Here the effective Boltzmann factor B{E) decays as an 

is a kinetic 



exponential of y/E. In particular, if £' = ^v? 
energy, one obtains distributions that exhibit exponen- 
tial tails in the velocity u |E|- This type of exponential 
behavior has been observed for stationary distributions 
of the complex Ginzburg Landau equation [s^], as well 
as in fusion plasma experiments [l9|. 



C. Stretched exponential tail 

As a generalization of the previous example, consider a 
density /(/?) which behaves as /(/?) ~ e^^^ , with c > 0, 
(5 < as /3 — > 0. This particular form of stretched ex- 
ponential implies a high-energy behavior of B{E) which 
also has the form of a stretched exponential. Indeed, 
solving the differential equation satisfied by Pe- 



[\nfm' = -cSP'-'=c\S\P'-'^E 



we find 



1/(5-1) 



(28) 



(29) 
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With this value oi Pe, the curvature of \nf{(3) is asymp- 
totically evaluated as follows: 



[ln/(/3ij)]" 



-c5iS - l)/3|r2 



= -c6{S-l) 



E 



(5-2)/(A--l) 



(30) 



Similarly, the Legendre transform of ln/(/5) is found to 
behave as 



^5/(5-1) 
(c|5|)l/(5-l) 

(c|^|)i/(5-i) 



(c|(5|)V(A-- 



(31) 



as -E — !■ cx). From Ea. (|15|l . we consequently obtain 

i?(i?)~i?(2-^)/(2^-2)ga£^/<^-i)^ (32) 

Stretched exponentials are relevant for observed dis- 
tributions in hydrodynamic turbulence 29], plasma ex- 
periments .30], as well as in granular media [sij. It is 
worth pointing out that the idea of superposing expo- 
nential factors to obtain a stretched exponential factor 
has been used previously by Palmer et al. |3R| to model 
anomalous relaxation dynamics; see also [33 for applica- 
tions of the same idea in the context of dissipative fluxes 
dynamics. 

D. Constant tail 

So far we have assumed that /(/?) ^ as /? 0. We 
next consider a case where /(/?) goes to some constant as 
/3 ^ 0. For this case, the differential equation defining 
Pe poses a problem because [ln/(/3)]' does not diverge 
as (3 approaches 0. To define (3e, we must then resort to 
find the largest value of $(/?, E) directly without the use 
of derivatives. As an example, consider the case where 
the /? ^ behavior of f{(3) is given by f{(3) ^ a, with 
a > 0. Then, 



Pe = argsup{-/?i? 4 
= argsup{— 



In a} 



(33) 



which implies that Pe — for all i? > 0. Consequently, 
B{E) must behave like a constant as i? ^ oo, since 



PEE + \iif{PE)^lna. 



(34) 



This is not a very interesting case physically because 
B{E) is not normalizable. Note that a constant asymp- 



totics for B{E) is also recovered if /(/?) 
c> 0. 



-cf3 



, a > 0, 



E. Log-normal distribution 

For our final example, we consider the case where P g 
(0, oo) is distributed according to the log-normal density 



/(/3) = -exp[-c(ln/3-6)2] 



(35) 



where a, h, and c are all positive constants. With this 
density, the integral defining the stationary distribution 
B{E) takes the form 



Jo P 



This is equivalent to 
B{E) =a 

J — ( 



^-c{y~-bf-Ee^ 



dy 



(36) 



(37) 



using the change of variables y = ln/3. At this point, 
we proceed as before to find the asymptotic behavior 
of B{E) as ^ oo by locating the saddle-point yE oi 
Ea. (|37fl which maximizes the exponent function 



<i>iy,E) = -c{y - hf - Eey 



(38) 



over all real values of y. The exact solution of ys can be 
found to be given by 



VE 



-W 



Ee'' 



(39) 



where W{x) is the Lambert or product-log function de- 
fined as the principal branch solution of the equation 
ye^ = X. The Laplace approximation of B{E) which 
results from this solution for j/b is 



liiB{E) 



'^{yE,E) 

'Ee'' 



-cW 



- Ee 



(40) 



A more manageable asymptotics can be found ana- 
lytically by expanding the exponential term in ^{y,E) 
around y = to first, second, or third order in y, and 
then find the maximum of the corresponding approxi- 
mation of ^{yE,E) to obtain an approximation of the 
Laplace asymptotics of B{E) found above. This method 
is described in Ref. 4Q], and yields surprising good re- 
sults already at third order in y despite the fact that 
yE — > — oo when E oo (recall that /? — > as oo 
and that y = ln/3). 



VI. INVERSE PROBLEM 

The inverse transform of the Laplace transform defin- 
ing B{E) has the form 



1 i'C+iOQ 

fiP) = 77- / B{E)e^''dE 



(41) 
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where c is an arbitrary constant lying in the domain of 
definition of B{E). From this integral, we are in a posi- 
tion to predict the /? — > or /3 ^ oo behavior that f{(3) 
needs to have in order for B{E) to behave according to 
some prescribed form. This is the inverse problem of the 
previous sections, namely: given a prescribed form for 
B{E), what is the behavior of /(/3)? 

We will not reflect much on this inverse problem be- 
cause the asymptotic methods that can be used to solve 
it are exactly the same as those described in the context 
of the direct problem. On the one hand, in the limit 
f} ^ 0, we may proceed just as in Sec. lIIII to expand the 
exponential term e'^^ in (|41l) to obtain a Taylor series for 

/(/?): 



/(/?) 



where 



— / B{E) [l + pE 



ao + ai/? + a2/3 -|- 



2m .L_i^ k\ 



dE 
(42) 

(43) 



The coefficient can be computed by rotating the com- 
plex integral shown above to the real line with the sub- 
stitution E — > lE. On the other hand, /3 — *■ oo asymp- 
totics for /(/3) can be found in just the same way as 
E oo asymptotics were found for B{E) using the 
Laplace method. At the level of the inverse Laplace 
transform integral defining /(/?), the application of this 
approximation method corrected with the Gaussian term 
yields 



fid) 



^l3E^+\nB{Ep) 



(44) 



for /3 ^ oo, where Efj is now found by solving the equa- 
tion 



[lnBiE)]'^-P 



(45) 



for E. This represents a valid approximation of /(/3) 
at large values of /3 (low-temperature limit) provided, 
as was the case for /(/?), that B{E) is unimodal and 
differentiable. It is important to notice that the energy 
value Efj taken as function of f3 is not the inverse function 
of In the direct problem, Pe is found for E oo, 
and in that limit we have seen that f3E — > 0. For the 
asymptotics of the inverse problem, however, we consider 
the limit /3 — > cx). 



VII. CONCLUSION 

We have analyzed the main asymptotic properties of 
general superstatistics which are convex superpositions 
of Boltzmann exponential factors. The saddle-point ap- 
proximation method turns out to be useful to treat this 
problem in full generality. In practice, our methods allow 
us to construct simple superstatistical models that may 
underlie an experimentally measured "fat tail" distribu- 
tion in a driven nonequilibrium system. In this context, 
we have seen that what universally determines the form 
of the high energy tail of observed distributions is the low 
f3 (high temperature) behavior of the mixing distribution 
/(/3) which defines at the bottom the superstatistical 
model. The spectrum of possibilities where our results 
can be applied is quite broad. It contains the power-law 
distributions of nonextensive statistical mechanics as a 
special case, but it is also relevant for stretched expo- 
nentials or lognormal superstatistics, as we have demon- 
strated. 
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